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Abstract Wc consider a basic model of the lossless interaction between a 
moving two-level atom and a standing-wave single-mode laser field. Classi- 
cal treatment of the translational atomic motion provides the semiclassical 
Hamilton-Schrodinger equations of motion which are a fivc-dimcnsional non- 
linear dynamical system with two integrals of motion. The atomic dynamics 
can be regular or chaotic (in the sense of exponential sensitivity to small 
variations in initial conditions and/or the system's control parameters) in 
dependence on values of the control parameters, the atom-field detuning and 
recoil frequency. We develop a semiclassical theory of the chaotic atomic 
transport in terms of a random walk of the atomic electric dipole moment 
u which is one of the components of a Bloch vector. Based on a jump-like 
behavior of this variable for atoms crossing nodes of the standing laser wave, 
we construct a stochastic map that specifies the center-of-mass motion. We 
find the relations between the detuning, recoil frequency and the atomic en- 
ergy, under which atoms may move in a rigid optical lattice in a chaotic 
way. We obtain the analytical conditions under which deterministic atomic 
transport has fractal properties and explain a hierarchical structure of the 
dynamical fractals. Quantum treatment of the atomic motion in a standing 
wave is studied in the dressed state picture where the atom moves in two 
optical potentials simultaneously. If the values of the detuning and a charac- 
teristic atomic frequency are of the same order, than there is a probability 
of nonadiabatic transitions of the atom upon crossing nodes of the standing 
wave. At the same condition exactly, we observe sudden changes (jumps) in 
the atomic dipole moment u when the atom crosses the nodes. Those jumps 
are accompanied by splitting of atomic wave packets at the nodes. Such a pro- 
liferation of wave packets at the nodes of a standing wave is a manifestation 
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of classical atomic chaotic transport. In particular, the effect of simultaneous 
trapping of an atom in a well of one of the optical potential and its flight in 
the other potential is a quantum analogue of a chaotic classical walking of 
an atom. At large values of the detuning, the quantum evolution is shown to 
be adiabatic in accordance with a regular character of the classical atomic 
motion. 



1 Short historical background 

The fundamental model for the interaction of a radiation with matter, com- 
prising a collection of two-level quantum systems coupled with a single-mode 
electromagnetic field, provides the basis for laser physics and describes a rich 
variety of nonlinear dynamical effects. The discovery that a single-mode laser, 
a symbol of coherence and stability, may exhibit deterministic instabilities 
and chaos is especially important since lasers provide nearly ideal systems 
to test general ideas in statistical physics. From the stand point of nonlinear 
dynamics, laser is an open dissipative system which transforms an external 
excitation into a coherent output in the presence of loss. In 1975 Haken [T] 
has shown that a single-mode, homogeneously broadened laser, operating on 
resonance with the gain center can be described in the rotating-wave ap- 
proximation by three real semiclassical Maxwell-Bloch equations which are 
isomorphic to the famous Lorenz equations. Some manifestations of a Lorenz- 
type strange attractor and dissipative chaos have been observed with different 
types of lasers. 

In the same time George Zaslavsky with co-workers [2] have studied inter- 
action of an ensemble of two-level atoms with their own radiation field in a 
perfect single-mode cavity without any losses and external excitations, which 
is known as the Dicke model [3J- They were able to demonstrate analyti- 
cally and numerically dynamical instabilities and chaos of Hamiltonian type 
in a semiclassical version of the Dicke model without rotating-wave approx- 
imation. It was the first paper that opened the door to study Hamiltonian 
atomic chaos in the rapidly growing fields of cavity quantum electrodynam- 
ics, quantum and atomic optics. Semiclassical equations of motion for this 
system may be reduced to Maxwell-Bloch equations for three real indepen- 
dent variables which, in difference from the laser theory, do not include losses 
and pump. Those equations are, in general, nonintegrable, but they become 
integrable immediately after adopting the rotating-wave approximation [3] 
that implies the existence of an additional integral of motion, conservation 
of the so-called number of excitations. Numerical experiments have shown 
that prominent chaos arises when the density of atoms is very large (approx- 
imately 10 20 cm 3 in the optical range [2]). The following progress in this field 
has been motivated, mainly, by a desire to find manifestations of Hamiltonian 
atomic chaos in the models more suitable for experimental implementations. 
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Twenty years after that pioneer paper, manifestations of Hamiltonian chaos 
have been found in experiments with kicked cold atoms in a modulated laser 
field. Nowdays, a few groups in the USA, Australia, New Zealand, Germany, 
France, England, Italy and in other countries can perform routine experi- 
ments on Hamiltonian chaos with cold atoms in optical lattices and traps 
(for a review see [12]). 

In this paper we review some results on theory of Hamiltonian chaos with 
a single two-level atom in a standing- wave laser field that have been obtained 
in our group in Vladivostok. In spite of we published with George only one 
paper on this subject [6], our work in this field has been mainly inspired by 
his paper [2] written in 1975 in Krasnoyarsk, Siberia. 



2 Introduction 

An atom placed in a laser standing wave is acted upon by two radiation forces, 
deterministic dipole and stochastic dissipative ones [7] . The mechanical action 
of light upon neutral atoms is at the heart of laser cooling, trapping, and 
Bosc-Einstein condensation. Numerous applications of the mechanical action 
of light include isotope separation, atomic lithography and epitaxy, atomic- 
beam deflection and splitting, manipulating translational and internal atomic 
states, measurement of atomic positions, and many others. Atoms and ions in 
an optical lattice, formed by a laser standing wave, are perspective objects for 
implementation of quantum information processing and quantum computing. 
Advances in cooling and trapping of atoms, tailoring optical potentials of 
a desired form and dimension (including one-dimensional optical lattices), 
controlling the level of dissipation and noise are now enabling the direct 
experiments with single atoms to study fundamental principles of quantum 
physics, quantum chaos, decoherence, and quantum-classical correspondence 
(for recent reviews on cold atoms in optical lattices see Ref. [5JI2])- 

Experimental study of quantum chaos has been carried out with ultracold 
atoms in <5-kickcd optical lattices [TUl EJ 02] . To suppress spontaneous emis- 
sion and provide a coherent quantum dynamics atoms in those experiments 
were detuned far from the optical resonance. Adiabatic elimination of the 
excited state amplitude leads to an effective Hamiltonian for the center-of- 
mass motion [13], whose 3/2 degrec-of- freedom classical analogue has a mixed 
phase space with regular islands embedded in a chaotic sea. De Brogile waves 
of <5-kicked ultracold atoms have been shown to demonstrate under appropri- 
ate conditions the effect of dynamical localization in momentum distributions 
which means the quantum suppression of chaotic diffusion [TUJ [TTJ [T^] . De- 
coherence due to spontaneous emission or noise tend to suppress this quan- 
tum effect and restore classical-like dynamics. Another important quantum 
chaotic phenomenon with cold atoms in far-detuned optical lattices is a chaos- 
assisted tunneling. In experiments |14| [5] ultracold atoms have been demon- 
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strated to oscillate coherently between two regular regions in mixed phase 
space even though the classical transport between these regions is forbidden 
by a constant of motion (other than energy). 

The transport of cold atoms in optical lattices has been observed to take 
the form of ballistic motion, oscillations in wells of the optical potential, 
Brownian motion [15] . anomalous diffusion and Levy nights [161 117] . The 
Levy flights have been found in the context of subrecoil laser cooling [16] in 
the distributions of escape times for ultracold atoms trapped in the potential 
wells with momentum states close to the dark state. In those experiments the 
variance and the mean time for atoms to leave the trap have been shown to 
be infinite. 

A new arena of quantum nonlinear dynamics with atoms in optical lattices 
is opened if we work near the optical resonance and take the dynamics of in- 
ternal atomic states into account. A single atom in a standing- wave laser field 
may be semiclassically treated as a nonlinear dynamical system with coupled 
internal (electronic) and external (mechanical) degrees of freedom [181 [T9l [20] . 
In the semiclassical and Hamiltonian limits (when one treats atoms as point- 
like particles and neglects spontaneous emission and other losses of energy), 
a number of nonlinear dynamical effects have been analytically and numer- 
ically demonstrated with this system: chaotic Rabi oscillations [HI HH1 HO] ; 
Hamiltonian chaotic atomic transport and dynamical fractals [21] [22] [23] [25] , 
Levy flights and anomalous diffusion [6] [20] [24] . These effects are caused by 
local instability of the CM motion in a laser field. A set of atomic trajectories 
under certain conditions becomes exponentially sensitive to small variations 
in initial quantum internal and classical external states or/and in the con- 
trol parameters, mainly, the atom-laser detuning. Hamiltonian evolution is a 
smooth process that is well described in a semiclassical approximation by the 
coupled Hamilton-Schrodingcr equations. A detailed theory of Hamiltonian 
chaotic transport of atoms in a laser standing wave has been developed in 
our recent paper [2"5] . 



3 Semiclassical dynamics 

3.1 Hamilton- Schro dinger equations of motion 

We consider a two-level atom with mass m a and transition frequency u a in a 
one-dimensional classical standing laser wave with the frequency uit and the 
wave vector kf. In the frame rotating with the frequency ujf, the Hamiltonian 
is the following: 

P 2 1 

H = 1 — h(uj a — u)f)a z — hfi (<x_ + <t_|_) coskfX. (1) 
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Here a± z are the Pauli operators which describe the transitions between 
lower, |1), and upper, |2), atomic states, fl is a maximal value of the Rabi 
frequency. The laser wave is assumed to be strong enough, so we can treat 
the field classically. Position X and momentum P operators will be consid- 
ered in section "Semiclassical dynamics" as c-numbers, X and P. The simple 
wavefunction for the electronic degree of freedom is 

|#(t)>=a(t)|2)+&(f)|l>, (2) 

where a and b are the complex- valued probability amplitudes to find the atom 
in the states |2) and |1), respectively Using the Hamiltonian ([T]), we get the 
Schrodingcr equation 

.da us a - u f 

i— = — -a — lib cos kfX, 

.db UJf -u a , 

i— = — o — J/a cos kfX. 

dt 2 1 

Let us introduce instead of the complex- valued probability amplitudes a and 
b the following real- valued variables: 

u = 2Rc(a&*), v = -2Im(a&*), z = \a\ 2 - \b\ 2 , (4) 

where u and v are a synchronized (with the laser field) and a quadrature 
components of the atomic electric dipole moment, respectively, and z is the 
atomic population inversion. 

In the process of emitting and absorbing photons, atoms not only change 
their internal electronic states but their external translational states change 
as well due to the photon recoil. In this section we will describe the transla- 
tional atomic motion classically. The position and momentum of a point-like 
atom satisfy classical Hamilton equations of motion. Full dynamics in the ab- 
sence of any losses is now governed by the Hamilton-Schrodinger equations 
for the real-valued atomic variables 

x = uj r p, t> ~ — it sin a;, u = Av, 
v = —Au + 2z cos x, z = —2v cos x, 

where x = kfX and p = P/hkf are normalized atomic ccnter-of-mass 
position and momentum, respectively. Dot denotes differentiation with re- 
spect to the dimensionless time r = fit. The normalized recoil frequency, 
u> T = hk 2 f/m a ri <C 1, and the atom-field detuning, A = (uif — w a )/fi, are 
the control parameters. The system has two integrals of motion, namely the 
total energy 

H= — v — ucosx z, 6 

2 2 
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and the Bloch vector u 2 + v 2 + z 2 = 1. The conservation of the Bloch vector 
length follows immediately from Eqs. 

Equations ([5]) constitute a nonlinear Hamiltonian autonomous system with 
two and half degrees of freedom which, owing to two integrals of motion, move 
on a three-dimensional hypersurface with a given energy value H. In general, 
motion in a three-dimensional phase space in characterized by a positive 
Lyapunov exponent A, a negative exponent equal in magnitude to the positive 
one, and zero exponent. The maximum Lyapunov exponent characterizes the 
mean rate of the exponential divergence of initially close trajectories and 
serves as a quantitative measure of dynamical chaos in the system. The result 
of computation of the maximum Lyapunov exponent in dependence on the 
detuning A and the initial atomic momentum p is shown in Fig. [TJ Color 
in the plot codes the value of the maximum Lyapunov exponent A. In white 




Fig. 1 Maximum Lyapunov exponent A vs atom-field detuning A and initial atomic mo- 
mentum pq: u> r = 10~ 5 , uo = zq = 0.7071, vq = 0. 

regions the values of A are almost zero, and the atomic motion is regular in 
the corresponding ranges of A and po- I n shadowed regions positive values of 
A imply unstable motion. 

Figure [T] demonstrates that the center-of-mass motion becomes unstable 
if the dimensionless momentum exceeds the value po ~ 300 that corresponds 
(with our normalization) to the atomic velocity v a ~ 3 m/s for an atom 
with m a ss 10~ 22 g in the field with the wavelength close to the transition 
wavelength A a ~ 800 nm. With these estimates for the atomic and lattice 
parameters and Q/2ir = 10 9 , one gets the normalized value of the recoil 
frequency equal tow r = 10~ J . The detuning A will be varied in a wide range, 
and the Bloch variables are restricted by the length of the Bloch vector. 
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3.2 Regimes of motion 

The case of exact resonance, A = 0, was considered in detail in Ref. |181 [24] . 
Now we briefly repeat the simple results for the sake of self-consistency. At 
zero detuning, the variable u becomes a constant, u = uq, and the fast (it, v, 
z) and slow (x, p) variables are separated allowing one to integrate exactly 
the reduced equations of motion. The total energy © is equal to Ho = H(u = 
u$,A = 0), and the atom moves in a simple cosine potential wocosx with 
three possible types of trajectories: oscillator-like motion in a potential well 
if Hq < uq (atoms are trapped by the standing- wave field) , motion along the 
separatrix if H = w , an d ballistic-like motion if H > u . The exact solution 
for the center-of-mass motion is easily found in terms of elliptic functions (see 



As to internal atomic evolution, it depends on the translational degree of 
freedom since the strength of the atom-field coupling depends on the position 
of atom in a periodic standing wave. At A = 0, it is easy to find the exact 
solutions of Eqs. J5]) 



where u = uq, and cos[x(r)] is a given function of the translational variables 
only which can be found with the help of the exact solution for x [HI [24] . 
The sign of v is equal to that for the initial value zq and xo is an integration 
constant. The internal energy of the atom, z, and its quadrature dipole- 
momcnt component v could be considered as frequency-modulated signals 
with the instant frequency 2cos[x(r)] and the modulation frequency oj r p(r), 
but it is correct only if the maximum value of the first frequency is much 
greater than the value of the second one, i. e., for \uj r p \ <C 2. 

The maximum Lyapunov exponent A depends both on the parameters w r 
and A, and on initial conditions of the system ([5]). It is naturally to expect 
that off the resonance atoms with comparatively small values of the initial 
momentum pq will be at once trapped in the first well of the optical potential, 
whereas those with large values of po will fly through. The question is what 
will happen with atoms, if their initial kinetic energy will be close to the 
maximum of the optical potential. Numerical experiments demonstrate that 
such atoms will wander in the optical lattice with alternating trappings in 
the wells of the optical potential and flights over its hills. The direction of the 
center-of-mass motion of wandering atoms may change in a chaotic way (in 
the sense of exponential sensitivity to small variations in initial conditions). 
A typical chaotically wandering atomic trajectory is shown in Fig. [2] 



us d]). 




(7) 
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Fig. 2 Typical atomic trajectory in the regime of chaotic transport: xq = 0, po = 300, 
20 = -1, uo = v = 0, ui r = 10~ 5 , A = -0.05. 



It follows from ([5]) that the translational motion of the atom at A ^ 
is described by the equation of a nonlinear physical pendulum with the 
frequency modulation 

x + Lo r u{r) sin a; = 0, (8) 
where u is a function of all the other dynamical variables. 



3.3 Stochastic map for chaotic atomic transport 

Chaotic atomic transport occurs even if the normalized detuning is very small, 
\A\ <C 1 (Fig. [T]). Under this condition, we will derive in this section approxi- 
mate equations for the center-of-mass motion. The atomic energy at |Z\| <C 1 
is given with a good accuracy by its resonant value H . Returning to the 
basic set of the equations of motion ([5]), we may neglect the first right-hand 
term in the fourth equation since it is very small as compared with the sec- 
ond one there. However, we cannot now exclude the third equation from the 
consideration. Using the solution ([7]) for v, we can transform this equation as 



u = ±A\fl—u 2 cosx, x = 2y cosxdr' + xo- (9) 

o 

Far from the nodes of the standing wave, Eq. ^ can be approximately inte- 
grated under the additional condition, \co r p\ <C 1, which is valid for the ranges 
of the parameters and the initial atomic momentum where chaotic transport 
occurs. Assuming cosx to be a slowly- varying function in comparison with 
the function cosx, we obtain far from the nodes the approximate solution for 
the w-component of the atomic dipole moment 

u » sin | ± — sinx + C) , (10) 
V 2 cos x J 
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where C is an integration constant. Therefore, the amplitude of oscillations 
of the quantity u for comparatively slow atoms (|w r p| <C 1) is small and of 
the order of |Z\| far from the nodes. 

At |Zi| = 0, the synchronized component of the atomic dipolc moment u is 
a constant whereas the other Bloch variables z and v oscillate in accordance 
with the solution ([7]). At \A\ ^ and far from the nodes, the variable u 
performs shallow oscillations for the natural frequency \A\ is small as com- 
pared with the Rabi frequency. However, the behavior of u is expected to 
be very special when an atom approaches to any node of the standing wave 
since near the node the oscillations of the atomic population inversion z slow 
down and the corresponding driving frequency becomes close to the resonance 
with the natural frequency. As a result, sudden "jumps" of the variable u are 
expected to occur near the nodes. This conjecture is supported by the nu- 
merical simulation. In Fig. [3] we show a typical behavior of the variable u 
for a comparatively slow and slightly detuned atom. The plot clearly demon- 




Fig. 3 Typical evolution of the atomic dipole-moment component u for a comparatively 
slow and slightly detuned atom: xq = 0, po = 550, vq = 0, uq = zo = 0.7071, u) r = 10 — 5 , 
A = -0.01. 



stratcs sudden "jumps" of u near the nodes of the standing wave and small 
oscillations between the nodes. 

Approximating the variable u between the nodes by constant values, we 
can construct a discrete mapping [23) 

u m = sin((9 sin </> rn + arcsinu rn _i), (11) 

where = \A\y/Tr/u) r p n ode will be called an angular amplitude of the jump, 
u m is a value of u just after the m-th node crossing, (j) m are random phases 
to be chosen in the range [0, 27r], and p no dc = \j2H/oj r is the value of the 
atomic momentum at the instant when the atom crosses a node (which is 
the same with a given value of the energy H for all the nodes). With given 
values of A, uj r , and p no dc, the map (|11[) has been shown numerically to 
give a satisfactory probabilistic distribution of magnitudes of changes in the 
variable u just after crossing the nodes. The stochastic map (fTTj) is valid 
under the assumptions of small detunings (|Z\| <^ 1) and comparatively slow 
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atoms (|w r p| «C 1). Furthermore, it is valid only for those ranges of the control 
parameters and initial conditions where the motion of the basic system 
is unstable. For example, in those ranges where all the Lyapunov exponents 
are zero, u becomes a quasi-periodic function and cannot be approximated 
by the map. 



3-4 Statistical properties of chaotic transport 

With given values of the control parameters and the energy if, the center- 
of-mass motion is determined by the values of u m (see Eq. ©). One can 
obtain from the expression for the energy ((SJ) the conditions under which 
atoms continue to move in the same direction after crossing a node or change 
the direction of motion not reaching the nearest antinode. Moreover, as in 
the resonance case, there exist atomic trajectories along which atoms move 
to antinodes with the velocity going asymptotically to zero. It is a kind of 
separatrix-like motion with an infinite time of reaching the stationary points. 

The conditions for different regimes of motion depend on whether the 
crossing number m is even or odd. Motion in the same direction occurs at 
(— l) rn+1 u m < H, separatrix-like motion — at (— l) m+1 u m = H, and turns — 
at (— l) m+1 u m > H. It is so because even values of m correspond to cos a; > 0, 
whereas odd values — to cos x < 0. The quantity u during the motion changes 
its values in a random-like manner (see Fig. [3]) taking the values which provide 
the atom either to prolong the motion in the same direction or to turn. 
Therefore, atoms may move chaotically in the optical lattice. The chaotic 
transport occurs if the atomic energy is in the range < if < 1. At if < 0, 
atoms cannot reach even the nearest node and oscillate in the first potential 
well in a regular manner (see Fig. [TJ. At if > 1, the values of u are always 
satisfy to the flight condition. Since the atomic energy is positive in the 
regime of chaotic transport, the corresponding conditions can be summarized 
as follows: at \u\ < if , atom always moves in the same direction, whereas at 
\u\ > if, atom either moves in the same direction, or turns depending on the 
sign of cos 2 in a given interval of motion. In particular, if the modulus of u 
is larger for a long time then the energy value, then the atom oscillates in a 
potential well crossing two times each of two neighbor nodes in the cycle. 

The conditions stated above allow to find a direct correspondence between 
chaotic atomic transport in the optical lattice and stochastic dynamics of the 
Bloch variable u. It follows from Eq. that the jump magnitude u m — u m —i 
just after crossing the m-th node depends nonlincarly on the previous value 
u m —i- For analyzing statistical properties of the chaotic atomic transport, it 
is more convenient to introduce the map for arcsin it m |23j 

9 m = arcsin u m = <9sin0 m + arcsin ii m _i, (12) 
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where the jump magnitude docs not depend on a current value of the variable. 
The map (|T2"j) visually looks as a random motion of the point along a circle 
of unit radius (Fig. @|. The vertical projection of this point is u m . The value 



Fig. 4 Graphic representation for the maps of u m and m = arcsin« m . H is a given 
value of the atomic energy. Atoms either oscillate in optical potential wells (trapping) or 
fly through the optical lattice (flight). 

of the energy H specifies four regions, two of which correspond to atomic 
oscillations in a well, and two other ones — to ballistic motion in the optical 
lattice. 

We will call "a flight" such an event when atom passes, at least, two 
successive antinodes (and three nodes). The continuous flight length L > 2it 
is a distance between two successive turning points at which the atom changes 
the sign of its velocity, and the discrete flight length is a number of nodes I 
the atom crossed. They are related in a simple way, L ~ irl, for sufficiently 
long flight. 

Center-of-mass oscillations in a well of the optical potential will be called 
"a trapping" . At extremely small values of the detuning, the jump magnitudes 
are small and the trapping occurs, largely, in the 27r-wide wells, i. e., in the 
space interval of the length 2tt. At intermediate values of the detuning, it 
occurs, largely, in the 7r-wide wells, i. e. in the space interval of the length it. 
Far from the resonance, \A\ > 1, trapping occurs only in the 7r-wide wells. 
Just like to the case of flights, the number of nodes /, atom crossed being 
trapped in a well, is a discrete measure of trapping. 

The PDFs for the flight Pr(1) and trapping Ptr(0 events were analytically 
derived to be exponential in a case of large jumps [23]. In a case of small 
jumps, the kind of the statistics depends on additional conditions imposed 
on the atomic and lattice parameters, and the distributions Pn(l) and Pt r (l) 
were analytically shown to be either practically exponential or functions with 
long power- law segments with the slope —1.5 but exponential "tails". The 
comparison of the PDFs computed with analytical formulas, the stochastic 
map, and the basic equations of motion has shown a good agreement in 
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different ranges of the atomic and lattice parameters [53]. We will use the 
results obtained to find the analytical conditions, under which the fractal 
properties of the chaotic atomic transport can be observed, and to explain 
the structure of the corresponding dynamical fractals. 

Since the period and amplitude of the optical potential and the atom-field 
detuning can be modified in a controlled way, the transport exponents of the 
flight and trapping distributions are not fixed but can be varied continuously, 
allowing to explore different regimes of the atomic transport. Our analytical 
and numerical results with the idealized system have shown that deterministic 
atomic transport in an optical lattice cannot be just classified as normal and 
anomalous one. We have found that the flight and trapping PDFs may have 
long algebraically decaying segments and a short exponential "tail" . It means 
that in some ranges of the atomic and lattice parameters numerical experi- 
ments reveal anomalous transport with Levy flights. The transport exponent 
equal to —1.5 means that the first, second, and the other statistical moments 
are infinite for a reasonably long time. The corresponding atomic trajectories 
computed for this time are self-similar and fractal. The total distance, that 
the atom travels for the time when the flight PDF decays algebraically, is 
dominated by a single flight. However, the asymptotic behavior is close to 
normal transport. In other ranges of the atomic and lattice parameters, the 
transport is practically normal both for short and long times. 



3. 5 Dynamical fractals 

Various fractal-like structures may arise in chaotic Hamiltonian systems [261 
[27] . In Ref. [22] [21] [24] [25] we have found numerically fractal properties of 
chaotic atomic transport in cavities and optical lattices. In this section we 
apply the analytical results of the theory of chaotic transport, developed in 
the preceding sections, to find the conditions under which the dynamical 
fractals may arise. 

We place atoms one by one at the point xq = with a fixed positive value of 
the momentum po and compute the time T when they cross one of the nodes 
at x = — 7r/2 or x = 3tt/2. In these numerical experiments we change the 
value of the atom-field detuning A only. All the initial conditions po = 200, 
zq = —1, uq = vq = and the recoil frequency tu r = 10~ 5 arc fixed. The exit 
time function T(A) in Fig. [5] demonstrates an intermittency of smooth curves 
and complicated structures that cannot be resolved in principle, no matter 
how large the magnification factor. The second and third panels in Fig. [5] 
demonstrate successive magnifications of the detuning intervals shown in the 
upper panel. Further magnifications reveal a self-similar fractal-like structure 
that is typical for Hamiltonian systems with chaotic scattering [26] [28] . The 
exit time T, corresponding to both the smooth and unresolved A intervals, 
increases with increasing the magnification factor. Theoretically, there exist 
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Fig. 5 Fractal-like dependence of the time of exit of atoms T from a small region in the 
optical lattice on the detuning A: po = 200, zo = — 1, Urj = i>o = 0. Magnifications of the 
detuning intervals are shown. 

atoms never crossing the border nodes at x = —ir/2 or x = 3tt/2 in spite 
of the fact that they have no obvious energy restrictions to do that. Tiny 
interplay between chaotic external and internal atomic dynamics prevents 
those atoms from leaving the small space region. 

Various kinds of atomic trajectories can be characterized by the number 
of times m atom crosses the central node at x = n/2 between the border 
nodes. There are also special separatrix-like trajectories along which atoms 
asymptotically reach the points with the maximum of the potential energy, 
having no more kinetic energy to overcome it. In difference from the separatrix 
motion in the resonant system (A = 0), a detuned atom can asymptotically 
reach one of the stationary points even if it was trapped for a while in a well. 
Such an asymptotic motion takes an infinite time, so the atom will never 
reach the border nodes. 

The smooth A intervals in the first-order structure (Fig. [5j upper panel) 
correspond to atoms which never change the direction of motion (m = 1) 
and reach the border node at x = 3tt/2. The singular points in the first-order 
structure with T = oo, which are located at the border between the smooth 
and unresolved A intervals, are generated by the asymptotic trajectories. 
Analogously, the smooth A intervals in the second-order structure (second 
panel in Fig. [5]) correspond to the 2-nd order (m = 2) trajectories, and so on. 

The set of all the values of the detunings, generating the separatrix-like 
trajectories, was shown to be a countable fractal in Refs. whereas the 

set of the values generating dynamically trapped atoms with m = oo seems 
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to be uncountable. The exit time T depends in a complicated way not only 
on the values of the control parameters but on initial conditions as well. 

In Fig. |5] [24| we presented a two-dimensional image of the time of exit T 
in the space of the initial atomic momentum pg and the atom-field detuning 
A. A self-similarity of this function is evident. 




Fig. 6 The scattering function in the regime of chaotic wandering. The time of exit T vs 
the detuning A and the initial momentum po- The function is shown in a shaded relief 
regime. 
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The length of all smooth segments in the m-th order structure in Fig. [5] 
is proportional to the number of atoms N(m) leaving the space [— tt/2, 3ir/2] 
after crossing the central node m times. An exponential scaling N(m) ~ 
exp(— 7m) has been found numerically with 7 ~ 1. The trapping PDFs, 
computed with the basic and reduced equations of motion at the detunings 
in the range shown in Fig. [5j have been found to have exponential tails. 
It is well known |26j that Hamiltonian systems with fully developed chaos 
demonstrate, as a rule, exponential decay laws, whereas the systems with a 
mixed phase space (containing islands of regular motion) usually have more 
slow algebraic decays due to the effect of stickiness of trajectories to the 
boundaries of such islands [27]. We have not found visible regular islands 
in our system at the values of the control parameters used to compute the 
fractal in Fig. [5] and we may conclude that the exponential scaling is a result 
of completely chaotic wandering of atoms in the space interval [— 7r/2, 3n/2] 
resembling chaotic motion in hyperbolic systems. 

The fractal-like structure with smooth and unresolved components may 
appear if atoms have an alternative either to turn back or to prolong the 
motion in the same direction just after crossing the node at x — tt/2. For 
the first-order structure in the upper panel in Fig. [5j it means that the in- 
ternal variable u of an atom, just after crossing the node for the first time 
(cos a; < 0), satisfies either to the condition u\ < H (atom moves in the same 
direction), or to the condition u\ > H (atom turns back). If u\ = H, then 
the exit time T is infinite. The jumps of the variable u after crossing the 
node are deterministic but sensitively dependent on the values of the control 
parameters and initial conditions. We have used this fact when introducing 
the stochastic map. Small variations in these values lead to oscillations of 
the quantity arcsinwi around the initial value arcsinuo with the angular am- 
plitude 0. If this amplitude is large enough, then the sign of the quantity 
ui — H alternates and we obtain alternating smooth (atoms reach the border 
x = 'Sir/ 2 without changing their direction of motion) and unresolved (atoms 
turns a number of times before exit) components of the fractal-like structure. 

If the values of the parameters admit large jump magnitudes of the variable 
it, then the dynamical fractal arises in the energy range < H < 1, i. e., at 
the same condition under which atoms move in the optical lattice in a chaotic 
way. In a case of small jump magnitudes, fractals may arise if the initial value 
of an atom u$ is close enough to the value of the energy H, i. e., the atom 
has a possibility to overcome the value u = H in a single jump. Therefore, 
the condition for appearing in the fractal T(A) the first-order structure with 
singularities is the following: 

I arcsinito — arcsini?| < 0. (13) 

The generation of the second-order structure is explained analogously. If 
an atom made a turn after crossing the node for the first time, then it will 
cross the node for the second time. After that, the atom either will turn or 
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cross the border node at x = —tt/2. What will happen depend on the value 
of it2- However, in difference from the case with m = 1, the condition for 
appearing an infinite exit time with m = 2 is u 2 = —H. Furthermore, the 
previous value u\ is not fixed (in difference from uq) but depends on the 
value of the detuning A. In any case we have u\ > H since the second-order 
structure consists of the trajectories of those atoms which turned after the 
first node crossing. In order for an atom would be able to turn after the second 
node crossing, the magnitude of its variable u should change sufficiently to 
be in the range u 2 < —H. The atoms, whose variables u could not "jump" so 
far, leave the space [— tt/2, 3tt/2]. The singularities are absent in the middle 
segment of the second-order structure shown in the second panel in Fig. [5] 
because all the corresponding atoms left the space after the second node 
crossing. The variable u 2 oscillates with varying A generating oscillations of 
the exit time. The condition for appearing singularities in the second-order 
structure is the following: 

2arcsiniJ < G. (14) 

With the values of the parameters taken in the simulation, we get the energy 
H = 0.2 + A/2. It is easy to obtain from the inequality (fT4| the approximate 
value of the detuning \A\ sa 0.0107 for which the second-order singularities 
may appear. In the lower panel in Fig.[5]one can sec this effect. No additional 
conditions are required for generating the structures of the third and the next 
orders. 

Inequality (|14|) is opposite to the inequality that determines the condi- 
tion for appearing power law decays in the flight PDF. Therefore, dynamical 
fractal may appear in those ranges of the control parameters where the Levy 
flights are impossible and vice versa. However, the trapping PDF may have 
a power law decay. Inequality (fT4|) in difference from (fl~3|) is strongly related 
with the chosen concrete scheme for computing exit times. It is not required 
with other schemes, say, with three antinodes between the border nodes. 



4 Quantum dynamics 

In this section we will treat atomic translational motion quantum mechan- 
ically, i. e., atom is supposed to be not a point particle but a wave packet. 
The corresponding Hamiltonian H has the form |T]) with X and P being the 
position and momentum operators. We will work in the momentum space 
with the state vector 

|#(*)>= / (a(P,t)\2)+b(pt)\l))\P)dP, (15) 



which satisfies to the Schrodinger equation 
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m^ = H\#). (i6) 

The normalized equations for the probability amplitudes have the form 

ia(p) = hcurp 2 - A)a{p) - \[b(p + 1) + 6(p- 1)], 

I I (17) 

ib(jp) = -{uj rP 2 + A)b(p) - -[a(p + 1) + o(p - 1)], 

with the same normalization and the control parameters as in the semiclas- 
sical theory. When deriving (|17|) . we used the following property of the mo- 
mentum operator P: 

C0S k f X\P) = \ (e ik f* + e~ ik f*) \P) = l(\P+ hk s ) + \P- hk f )) . (18) 

Equations (|17[) are an infinite-dimensional set of ordinary differential complex- 
valued equations of the first order with coupled amplitudes a{j> ± n) and 
b(p± m). To characterize the internal atomic state, let us introduce the fol- 
lowing variables; 



u(r) = 2 Rc / dx [a(x, r)b* (x, r)] , 
^)- 2] m/* W ,,r^r)], (19) 
z{t)s ! 4t[|o(i,t)| 2 -|6(i,t)| 2 ], 



which are quantum versions of the Bloch components (jU, and we denote 
them by the same letters. 



5 Dressed states picture and nonadiabatic transitions 

Interpretation of the atomic wave-packet motion in a standing-wave field is 
greatly facilitated in the basis of atomic dressed states which are eigenstates 
of a two-level atom in a laser field. The adiabatic dressed states 

\ + ) A =sin0|2) +cos0|l), \-) A =cos<9|2) -sin0|l), 



A If A \ 2 (20) 

tan6> = \ + 1 

2 cos x V V 2 cos x ) 



are eigenstates at a nonzero detuning. The corresponding values of the 
quasienergy are 
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Fig. 7 Resonant Eq ' and nonresonant ' potentials for an atom in a standing wave. 
The optical Stern-Gcrlach effect in the resonant potential is shown: splitting of an atomic 
wave packet launched at the node of the wave (x$ = 7r/2, po = 0). The wave packet, placed 
initially at the antinode (xq = 0, po = 0), appears to be simultaneously at the top of E^~^ 
and the bottom of Eq ' potentials. Its | + )-component slides down both the sides of Eg + ^ 
and the | — )-component oscillates at the bottom of Eq ' . 




(21) 

Figure [Jj shows a spatial variation of the quasienergies along the standing- 
wave axis. It follows from Eqs. (f2"0)) and (|21l) that, in general case, atom moves 

in the two potentials E^' simultaneously. 

At exact resonance, A = 0, the dressed states have the simple form 

|+> = ^(|1> + |2», |-) = -L(|l)-|2» (22) 

and are called diabatic states. The resonant potentials, E^ = ±cosx, cross 
each other at the nodes of the standing wave, x = n/2+irm, (m = 0, ±1, . . .). 
What will happen if we place the centroid of an atomic wave packet exactly 
at the node, Xq = 7r/2, in the ground state |1) and suppose its initial mean 
momentum to be zero, po = 0? The initial ground state is the superposition 
of the diabatic states: |l) = (|+) + | — ))/v / 2- One part of the initial wave 
packet at the top of the potential Eq will start to move to the right under 
the action of the gradient force = —dE^/dx = sinx, and another 

one — to the left to be forced by F^ = — sinx (see Fig. [7]). It is the well- 
known optical Stern- Gerlach effect HHUZHSD]- If the maximal expected value 
of the atomic kinetic energy does not exceed the potential one, the atom will 
be trapped in the potential well. Two splitted components of the initial wave 
packet will oscillate in the well with the period of oscillations 
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The wave packet, with p a = 0, placed at the antinode, say, at xo = 0, is 
simultaneously at the top of the potential Eq and at the bottom of Eq '. 
Therefore, its | + ) -component will slide down the both sides of the potential 
curve Eq + \ and the | — )-component will oscillate around x — (see Fig.[7|). 

Out off resonance, A ^ 0, the atomic wave packet moves in the bipotential 
E A ±] (HI]). The distance between the quasienergy curves is minimal at the 
nodes of the standing wave and equal to A (see Fig. [7]) . The spatial period 
and the modulation depth of the resonant potentials Eq^ are twice as much 
as those for the nonresonant potentials E A . 

The probability of nonadiabatic transitions between the dressed states 
| + ) A and | — ) A can be estimated in a simple way. The time of flight over 
a short distance Sx in neighbourhood of a node is 6x/tjj r p no dc- If the time of 
transition between the quasienergy levels, 2 /A, is of the order of the flight 
time, the transition probability is close to 1. It is easy to get the characteristic 
frequency of atomic motion from that condition [7] 

^0 = \M-Pnodc, (24) 

where p no dc is a value of the momentum in the vicinity of a node. 

Depending on the relation between A and Aq, there are three typical cases. 

1. If \A\ Aq, the nonadiabatic transition probability between the states 
| + ) A and | — } A upon crossing any node is close to 1. However, the dia- 
batic states | + ) and | — ) arc not mixed, and atom moves in one of optical 
resonant potentials. 

2. If \A\ ~ Aq, the atom may or may not undergo a transition upon crossing 
any node from one of the nonresonant potentials to another one with the 
probabilities of the same order. 

3. If \A\ ^> Aq, the nonadiabatic transition probability is exponentially small, 
and atom moves in one of the nonresonant potentials. 



5.1 Wave packet motion in the momentum space 

The atom at r = is supposed to be prepared as a Gaussian wave packet in 
the momentum space 



ao(p)=0, bo(p) = 



1 



\J^A V 



exp 



2{A P ) 2 



i(p~Po)x 



(25) 



with the momentum width Ap —10 corresponding to the spatial width AX = 
A//407T that is much smaller than the optical wavelength Xf. We compute 
the probability to find a two-level atom at the moment of time r with the 
momentum p 
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W(p,T) = \a(p,T)\ 2 + \b(p,T)f 



(26) 



by integrating Eqs.([T71) with the initial condition (|2"5)l . The recoil frequency, 
oj r = 10 ~ 5 , is fixed and the centroid of the wave packet is placed at the 
antinodc xq = 0, in all the numerical experiments. 



5.1.1 Adiabatic evolution at exact resonance 




Fig. 8 Time dependence of the momentum probability function W(p, r) for a ballistic 
atom at resonance prepared initially in the ground state (A = 0, u) r = 10~ 5 , xq = 0, 
Po = 800). 



At exact resonance, A = 0, the wave functions of the diabatic states 
| + ) and | — ) evolve independently, each one evolves in its own potential 
£ , q + ' and Eq ^ , respectively. The atom, prepared initially in the ground state 
|1) = (| + ) + |-))/v / 2with the mean initial momentum po = 800, will start 
to move from the top of Eq + ^ and the bottom of Eq ' potentials (see FigJTJ) . 
Thus, the initial wave packet will split into two components | + ) and | — ). 
Time evolution of the probability function (|26]l for each of the components 
is shown in FiglHl Pay, please, attention that the values of p on this and 
similar plots increase downwards. Color in this figure codes the values of 
W(p, t). The | + )-component (the lower trajectory in the figure) slides down 
the curve E^ and, therefore, moves with an increasing velocity up to the 
next antinodc at x = n, and then it slows down approaching the antinodc 
at x = 2tt. The atom moves in the positive direction of the axis x and the 
process repeats periodically with the period Tq + ' = 2-K/oj r p\ ^ ~ 690, where 

Po~2-k i s a mean momentum of the | + )-component upon the atomic motion 
between and 2ir. 
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The | — )-component (the upper trajectory in FigjS]) moves upward the 
potential curve Eq ^ and slows down up to reaching the top of Eq ' at x = ir. 
Then it moves with an increasing momentum up to x = 2ix. Since the mean 
momentum of the | — ) -component is smaller than that of the | + ) one, the 
corresponding period is longer, Tq ~ 980. 

5.1.2 Proliferation of wave packets at the nodes of the standing 
wave 



440- 
P 




Fig. 9 Proliferation of atomic wave packets at the nodes of the standing wave at the de- 
tuning A = 0.05. The atom is prepared initially in the dressed state | + }. Other conditions 
are the same as in Fig[8] 

New features in propagation of atomic wave packets through the standing 
wave appear under the condition A ~ Aq. Using the semiclassical expres- 
sion for the total atomic energy ©, let us estimate the value of the atomic 
momentum at the nodes of the standing wave if the detuning is not large, 
\A\ 1. If the atom is prepared initially in the state | +), i.e., uq = 1, 
zq = 0, and xq = then we have H = Hq = 2.2 at po = 800. Since the total 
energy is a constant, we get immediately from Eq. © 

p no do =i y/2H/uj r ~ 665. (27) 

Using the same formula ([5]), we get the values of the minimal and maximal 
momenta if the atom starts to move with the initial mean momentum po = 
800: Pmin ~ ^/2(H Q - l)/u r 490 and p max ~ ^2{H + l)/w r ~ 800. 

The formula (|2^|) gives us the value of the characteristic frequency under 
the chosen conditions, Aq ~ 0.08. We fix A = 0.05 in this section, so A ~ Aq. 



22 



S.V. Prants 



The initial state | + ) is the following superposition of the adiabatic states: 

I +) = ^=[(cos0 + sin 0)| + ) A + (cos0 - sin 0)| -) J. (28) 

With the help of (|2~T|) we can estimate the mixing angle at A = 0.05 to 
be equal to 6 ~ —n/4. Then it follows from ((28} that | + ) ~ | — ) A , i. e., 
practically all the wave packet is initially at the bottom of the potential E A 
(Fig. [7j). Figure [9] demonstrates that the wave packet really slows down, and 
its centroid intersects the node x = n/2 at t± ~ 215. Under the condition 
A ~ Aq, the atom has a probability to change the potential for another 
one upon crossing a node and a probability to stay in its present potential. 
This is exactly what we see in fig. [9j the wave packet splits at the node 
x = 7r/2 with the | — )-component climbing over the potential E^ (see the 
upper trajectory in this figure) and the | + )-component sliding down the 
curve E { A } with an increasing momentum (see the lower trajectory). Just 
after crossing the node, the most part of the probability density moves in 
the potential E A because the corresponding probability is larger. The | + )- 
component increases its velocity upon approaching the antinode at x = n 
and then slows down up to the second node at x = 3tt/2 where it splits 
into two components at ~ 640. After that, one of the components will 
move in the potential E A decreasing the velocity up to the next antinode 
at x = 2ir, and the other one will move in E A ' increasing its velocity in the 
same space interval. The probability density of this | — )-component is only 
a few percents, and we draw a solid curve along this trajectory in order to 
visualize the motion. 
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Fig. 10 The same as in Fig[9]but for the atom prepared initially in the ground state. 
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The | — ) -component of the packet, splitted after crossing the first node 
at x — 7r/2, has a smaller mean momentum than the | + )-one. Therefore, it 

reaches the second node later, at — 800, where it splits into two parts: 
the upper | + )-component will move in the potential and the lower 

| — )-one — in E A ~ . Such a proliferation of atomic wave packets takes places 
upon crossing all the next nodes of the standing wave. 

The moment of time Tn , when the centroids of the | ± )-components cross 
the n-th node, can be estimated by the simple formula (we suppose that the 
centroid of the atomic wave packet was at x = at r = 0): 

"rP { n % n T 7 \ ±] = (2« - 1) \ , n = 2, 3, . . . . (29) 

where p„_ i n is a mean momentum of the | ± ) -components upon their 
movement between (n — l)-th and n-th nodes. This quantity for the | — )- 
component, moving between x = and x = 7r/2, is Pq ^ = (po +p n odc)/2 — 
732.5. So, the centroid of this wave packet crosses the first node at t[~^ ~ 214. 
The lower | + )-component crosses the second node at x = 3ir/2 at t% — 
642. For the upper | — )-component we get p\ ^ = (p ri odc +Pmin)/2 — 577.5 

and T2 ~ 815. All the other moments of time, t« , can be estimated in 
the same way. The estimates obtained fit well the numerical data (see Fig^J) ■ 
The interference fringes on the upper trajectory at r ~ 1000 and p ~ 500 
and on the lower one at r ~ 900 and p ~ 800 reflect the fine-scale splitting 
of the corresponding wave packets. 

Let us now compute the probability map for the atom prepared initially 
in the ground state |1) which has the following form in the adiabatic state 
basis: 

|1) =cos0| + ) A -sin0|-) A , (30) 

It follows from ([2Tj) that (|30|) is almost a 50%~50% superposition of the 
I + ) A and I — ) A states. All the other conditions are assumed to be the same 
as before. The atomic wave packet splits from the beginning into two com- 
ponents with the I + )-one sliding down the curve E A (the lower trajectory 
in Fig. [TO]) and the | — )-one climbing over the potential E A ^ (the upper 
trajectory). Each of the components splits at the first node with a small 
time difference between the events. The subsequent proliferation of the wave 
packets occurs for the upper and lower parts of the probability density inde- 
pendently on each other in accordance with the same scenario as described 
above. In difference from the preceding case, the atom, prepared initially in 
the ground state, acquired the values of the momentum that are larger then 
the initial momentum po = 800. 

The nonadiabatic transitions are accompanied by drastic changes in the in- 
ternal state of the atom which is characterized by the values of the synphased 
component of the electric dipole moment u and the population inversion z. 
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Fig. 11 Time dependence of the dipole moment u and the population inversion z at the 
same conditions as in Fig. [9] 

In Fig. [IT] we demonstrate their behavior for the atom prepared initially in 
the state | + ) . Both the variables change their values abruptly in the time 
intervals with the centers at r ~ 215, 640 and 815, i. e., when the centroids 
of the atomic wave packets cross the first two nodes. 

5.1.3 Adiabatic motion at large detunings 

For comparison with the results of the preceding section, we demonstrate in 
Fig. [12] the evolution of the momentum distribution function W(p, r) with 
the ground initial state at A = 2 and the other same conditions as in the 
preceding section. The detuning A = 2 is large as compared to the character- 
istic frequency Ao ~ 0.09 that is estimated from (|24|) at po = 800. It follows 
from (|2T)]) and (|2"Tj) that at A — 2 the initial state |1) is a superposition of 
approximately 70% of the state | + }* and ~ 30% of the state | — So the 
main part of the initial packet begins to move in the potential E A increas- 
ing the momentum upon approaching the node at x — n/2, and the other 
part moves in E A ^ decreasing the momentum in the same space interval (see 
Fig. ri"2"T) . Upon crossing the nodes, the probability of transition between the 
states | ± ) A is small if \A\ Aq, and each of the component will continue 
to move in its own potential. The process is repeated and we see the periodic 
variations of the mean momentum of each of the components. The same pic- 
ture is observed if we take the state | + ) = (|1) + |2))/\/2 as the initial one. 
At A = 2, the state | + ) is a mix of 70% of | - ) A and 30% of | + ) A , so the 

main part of the initial | 4- ) wave packet will move in the potential E A ^ . The 
evolution of the internal atomic variables z and u is shown in Fig. 1131 There 
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are no jumps of z and u when the atom crosses nodes. Instead of that, we see 
fast oscillations of those variables when the atom crosses the first antinodcs. 




400 800 1200 1600 2000 



Fig. 12 Adiabatic evolution of the momentum probability function W(p, t) for a ballistic 
atom at the large detuning A = 2. 




Thus, at \A\ 3> Z\rj, there are no nonadiabatic transitions due to the cor- 
responding small probability and, therefore, no proliferation of wave packets 
at the nodes. The evolution of the atomic wave packet is adiabatic. 
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5.1.4 An atom can fly and be trapped simultaneously 

An intriguing effect of simultaneous trapping of an atom in a well of the 
optical potential and its ballistic flight through the optical lattice is observed 
at comparatively small values of the detuning. Let us prepare an atom in 
the ground state |1) with such a mean initial value of the momentum po 
that its | — ) -component would not be able to overcome the barrier of the 
potential E A ' but its | + ) -component would have a sufficient kinetic energy 
to overcome the barrier of the potential. Now one could expect periodic 
oscillations in the first well of the potential E^ A ^ and a simultaneous ballistic 

flight in the E^ potential with a proliferation of wave packets of the | + )- 
component at the nodes of the standing wave. 

Figure [141 demonstrates this effect at pq = 300, A = —0.05 and the same 
other conditions as before. We see that the momentum of the | — )-component 
(the upper trajectory in this figure) oscillates in the range (300, —300), and 
this component is trapped in the first well (— tt/2 < x < tt/2). Whereas 
the | +) -component moves in the positive direction splitting at each node. 
Estimates of the period of oscillations of the | — )-component, T ~ 2240, with 
the help of ([2U]) and of the time when the centroid of the | + ) -component 
crosses the first node, t[ ~ 380 (formula ([29]) ). fit well the data in Fig. H4l 
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Fig. 14 Effect of simultaneous trapping of an atom in a well of the optical potential and 
its flight through the wave. The ground initial state, A = —0.05, po = 300. 
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6 Quantum-classical correspondence and manifestations 
of dynamical chaos in wave-packet atomic motion 

Dynamical chaos in classical systems is characterized by exponentially fast 
divergence of initially close trajectories in a bounded phase space. Such a 
behavior is possible because of the continuity of the classical phase space 
whose points (therefore, classical system's states) can be arbitrary close to 
each other. The trajectory concept is absent in quantum mechanics whose 
phase space is not continuous due to the Heiscnberg uncertainty relation. 
The evolution of an isolated quantum system is unitary, and there can be no 
chaos in the sense of exponential sensitivity of its states to small variations 
in initial conditions. What is usually understand under "quantum chaos" is 
special features of the unitary evolution of a quantum system in the range of 
its parameter values and initial conditions at which its classical analogue is 
chaotic. 

The question "what happens to classical motion in the quantum world" 
is a core of the problem of quantum-classical correspondence. In spite of 
years of discussions from the beginning of the quantum era, it is still unclear 
how classical features appear from the underlying quantum equations. It is 
especially difficult to specify what happens to classical dynamical chaos in 
the quantum world [31j E2 El El E3 EH [37] . The interest to the problem 
of "quantum chaos" is motivated by our desire to understand the quantum 
origin of the observed classical chaos. 

In this section we establish a correspondence between the quantized mo- 
tion of a two-level atom in a standing laser wave and its semiclassical ana- 
logue considered in the third section. Semiclassical equations ^ represent a 
nonlinear dynamical system with positive values of the maximal Lyapunov 
exponent in a wide range of the initial conditions and control parameters to r 
and A. In other words, trajectories in the five-dimensional phase space are 
exponentially sensitive to small variations in initial conditions and/or param- 
eters in those ranges. That local dynamical instability is a reason for chaotic 
Rabi oscillations and chaotic motion of the atomic center of mass discussed 
in the third section. In particular, it has been found that an atom is able 
to walk chaotically in a strictly periodic optical lattice without any noise or 
other random processes (see Fig. The chaotic behavior is caused by jumps 
of the electric-dipole moment u at the nodes of the standing wave (Fig. [3]) . 
It follows from Eqs. ([5]) that this quantity governs the atomic momentum. A 
stochastic map for the quantity u (|TTj) allowed to derive analytic expressions 
for probability density functions of the atomic trapping and flight events that 
have been shown to fit well numerical simulation |23j . 

It has been shown that sudden changes in the behavior of u take place when 
we quantized the atomic motion (see Fig. ITTj) under the condition A ~ Aq. 
Those changes are more smooth than the jumps of u in the semiclassical case 
because a delocalized wave packet crosses a node for a finite time interval. 
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The quantum analysis provides a clear reason for those jumps at A ~ A , 
namely, it is nonadiabatic transitions between the quasienergy states | + ) A 
and | — ) A which occur when an atom crosses any node of the standing wave. 
Those jumps are accompanied by splitting of wave packets at the nodes. 
We may conclude that the proliferation of wave packets at the nodes of the 
standing wave is a manifestation of classical chaotic transport of an atom 
in an optical lattice that has been shown in Refs. [2TJ [24j [23] to take place 
in exactly the same ranges of initial conditions and control parameters. In 
particular, the effect of simultaneous trapping of an atom in a well of the 
optical potential and its flight in the same potential (Fig. I14[) is a quantum 
analogue of a chaotic walking of an atom shown in Fig. [2] 

In conclusion we would like to discuss briefly the role of dissipation. We 
did not take into account any losses in our treatment. Coherent evolution of 
the atomic state in a near-resonant standing- wave laser field is interrupted by 
spontaneous emission events at random moments of times. The semiclassical 
Hamiltonian evolution between these events has been shown to be regular or 
chaotic depending on the values of the detuning A and the initial momen- 
tum po. We stress that dynamical chaos may happen without any noise and 
any modulation of the lattice parameters. It is a specific kind of dynamical 
instability in the fundamental interaction between the matter and radiation. 

Dissipative transport of spontaneously emitting atoms in a ID standing- 
wave laser field has been studied in detail in Ref. [38] in the regimes where 
the underlying semiclassical Hamiltonian dynamics is regular and chaotic. A 
Monte Carlo stochastic wavefunction method was applied to simulate semi- 
classically the atomic dynamics with coupled internal and translational de- 
grees of freedom. It has been shown in numerical experiments and confirmed 
analytically that chaotic atomic transport can take the form either of bal- 
listic motion or a random walking with specific statistical properties. The 
character of spatial and momentum diffusion in the ballistic atomic trans- 
port was shown to change abruptly in the atom-laser detuning regime where 
the Hamiltonian dynamics is irregular in the sense of dynamical chaos. A 
clear correlation between the behavior of the momentum diffusion coefficient 
and Hamiltonian chaos probability has been found. 

What one could expect if spontaneous emission would be taken into consid- 
eration with our fully quantum equations of motion? Any act of spontaneous 
emission interrupts a coherent evolution of an atom at a random time mo- 
ment and is accompanied by a momentum recoil and a sudden transition of 
the atom into the ground state which is a superposition of the dressed states. 
The coherent evolution starts again after that. A collapse of the atomic wave 
function and a splitting of atomic wave packets are expected just after any 
spontaneous emission event. That additional splitting of wave packets at ran- 
dom time moments, besides of their proliferation at the nodes of a standing 
wave at A ~ Aq, can improve the quantum-classical correspondence in the 
regime of Hamiltonian chaos. 
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